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ABSTRACT 

We take advantage of the sensitivity and resolution of the Herschel Space Observatory at 100 and 
160 /im to directly image the thermal dust emission and investigate the infrared luminosities (£m) 
and dust obscuration of typical star-forming (L*) galaxies at high redshift. Our sample consists of 
146 UV-selected galaxies with spectroscopic redshifts 1.5 < z S pcc < 2.6 in the GOODS-North field. 
Supplemented with deep Very Large Array (VLA) and Spitzer imaging, we construct median stacks 
at the positions of these galaxies at 24, 100, and 160 /im, and 1.4 GHz. The comparison between 
these stacked fluxes and a variety of dust templates and calibrations implies that typical star-forming 
galaxies with UV luminosities Ltjv <] 10 10 L at z ~ 2 are luminous infrared galaxies (LIRGs) with 
a median L\r = (2.2 ± 0.3) x 10 u Lq. Their median ratio of Lir to rest-frame 8 fim luminosity 
(ig) is I/ir,/Z/8 = 8.9 ±1.3 and is « 80% larger than that found for most star-forming galaxies at 
z < 2. This apparent redshift evolution in the Ltr/Ls ratio may be tied to the trend of larger infrared 
luminosity surface density for z > 2 galaxies relative to those at lower redshift. Typical galaxies at 
1.5 < z < 2.6 have a median dust obscuration Lir/Ltjv = 7.1 ± 1.1, which corresponds to a dust 
correction factor, required to recover the bolomctric star formation rate (SFR) from the unobscured 
UV SFR, of 5.2 ± 0.6. This result is similar to that inferred from previous investigations of the UV, 
Ha, 24 /im, radio, and X-ray properties of the same galaxies studied here. Stacking in bins of UV 
slope (/3) implies that L* galaxies with redder spectral slopes are also dustier, and that the correlation 
between /3 and dustiness is similar to that found for local starburst galaxies. Hence, the rest-frame 
~ 30 and 50 /jm fluxes validate on average the use of the local UV attenuation curve to recover the 
dust attenuation of typical star-forming galaxies at high redshift. In the simplest interpretation, the 
agreement between the local and high redshift UV attenuation curves suggests a similarity in the dust 
production and stellar and dust geometries of starburst galaxies over the last 10 billion years. 
Subject headings: galaxies: high-redshift — infrared: galaxies — ISM: dust, extinction 
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1. INTRODUCTION 

A key aspect of measuring total star formation rates is 
to assess how a galaxy's spectral energy distribution is 
modulated by the effects of interstellar dust. The ability 
to quantify dust attenuation at high redshift is hindered 
by the limited dynamic range in luminosity at wave- 
lengths that are sensitiv e to dust emission While the 
Lyman Break technique (jSteidel et al.lfl995h has proven 
to be the most effective method of identifying galaxies 
across a large dynamic range in luminosity and look- 
back time, it requires observations in the rest-frame UV: 
the massive stars giving rise to the UV continuum also 
produce much of the dust that attenuates this contin- 
uum. The diminution of UV light is compounded by the 
fact that much of this radiation is re-emitted in the far- 
infrared where current instrumental sensitivity is insuf- 
ficient to directly detect typical star-forming galaxies at 
z > 2. It has therefore become common to use local rela- 
tions between monochromatic and bolometric luminosity 
to infer the star formation rates and dust attenuation of 
high-redshift galaxies. 

Prior to the e ra of large-scale multi- w a velength surveys 
like GOODS (jDickinson et al.l 120031 : iGiavalisco etall 
2004), inferring the dust attenuation at high redshift 
commonly entailed using relations between the variation 
of the UV continuum sl ope with dustiness as found for lo- 
cal starburst galaxies (jMeurer et al.lll999t ICalzetti et al.l 



l2000t iHeckman et al.l 119981 ) . However, such correlations 
were previously untested at high redshift. The correla- 
tion between the redness of the UV slope (/3) and dust 
attenuation has limited applicability when examined over 
a larger range in galaxy stellar population and luminos- 
ity. Galaxies with older and less massive stars contribut- 
ing significantly to th e UV em i ssion can also exhibit a 
redder spectral slope (jCalzettil Il997t Kong et all 12004 
iBuat et aLll2005f) . Further, the UV slope decouples from 
extinction for very luminous galaxies where virtually all 
of the star formation is obscured, such as is the case for 
low redshift ultra-lum inous infrared galaxies (ULIRGs; 
iGoldader et al.l 12002ft . Nonetheless, the local trend be- 
tween UV slope and dustiness is still applied widely to 
galaxies at high redshift; this correlation is often the only 
means by which one can infer the dust attenuation of 
galaxies at z > 3, where the dust emission from a typi- 
cal galaxy may be several orders of magnitude below the 
sensitivity limits of current infrared instruments. 

Incremental progress in determining dust attenua- 
tion at high redshift was achieved wit h the first ultra- 
deep radio and X-ra y surveys (e.g., I Richards] 120001 : 
lAlexander et al.1 12003ft . Though the sensitivity (even 
at GOODS-depth) at these wavelengths was insuffi- 
cient to individually detect typical star-forming galax- 
ies at z > 2, these surveys did allow for estimates 
of the mean dust attenuation based on stacking anal- 
yses. Initial studies suggested a general agreement be- 
tween UV and radio/ X-ray inferences o f dust attenu- 
ation at z ~ 2 — 3 (|Nandra et al.l |2002|: ISeibert et all 
2002t iReddy fe Steidell 12004 IReddv et all l2005l~[2006b; 
Daddi et al.l l2007t iPannella et alll2009ft . 



The launch of the Spitzer Space Telescope enabled 
the first direct detection of the dust emission in non- 
lensed L* galaxies at z ~ 2; at these redshifts, 
Spitzer's MIPS 24 /zm band is sensitive to the strongest 
dust emission feature (at 7.7 (im) in star-forming galax- 
ies. This feature arises from the stochastic UV photo- 
heating of small du s t grains and hydro carbons (e.g., 
iPuget fe Leeerl 119891: iTielens et all 11999ft and is found 
to correlate with the UV rad i ation from OB stars (e.g., 
iForster Schreiber et all 12004 iRoussel et all l2001ft. al- 
beit with si gnificant scatter (e.g.. iKennicutt et alll2009l: 
Hogg et all 120051 iHelou et al.l 1200 fl lEngclbracht et all 
20051: iNormand et al.lll995ft . Several studies have demon- 



strated that the dust attenuation inferred from the UV 
spectral slope, j3, for typical z ~ 2 galaxies is in gen- 
eral agreement with that inferred from MIPS 2 4 /xm (e.g., 
IReddv et all l2006bl l2010l: iDaddi et al.l 12007ft . Regard- 
less, local studies of resolved galaxies have emphasized 
the complexity of the 8 (im emission and the need to 
combine it with other obscured (IR) and unobscured 
(UV, Ha) tracers of star formation in order to obtain 
the most reliable meas ure of the bolometric luminosity 
(jKennicutt et al.ll2009ft . 

To this end, we take advantage of the improved 
sensitivity of Hers c hel/PACS at 100 and 160 /im 
(jPilbratt et all 120101 : iPoglitsch et~aT1 12010ft to measure 
directly for the first time the thermal dust emission from 
a large sample of typical (L*) star-forming galaxies at 
z ~ 2. We make use of the deep PACS data made pos- 
sible by the GOODS- Herschel Open Time Key project 
(PI: D. Elbaz). We supplement these data with deep 
Very Large Array (VLA) 1.4 GHz continuum imaging 



in the GOODS-North field (jMorrison et all 12010ft . The 
primary aim is to constrain the average infrared lumi- 
nosities of galaxies at high redshift, particularly those 
selected by their rest-frame UV colors, and to compare 
them to UV-based inferences. We begin in Section [2] by 
discussing the UV-selected galaxies and previous efforts 
to measure their stellar populations and dust attenua- 
tion. We also provide a brief description of the Herschel 
and radio data. Our stacking method and stacking sim- 
ulations are described in Section [31 followed by a discus- 
sion of the dust spectral energy distribution (SED) fits 
to the stacked fluxes and the total infrared luminosities 
in Section |4] In Section [5] we proceed to compare the 
Herschel-b&sed inferences of the dust attenuation with 
that inferred from the UV slope and discuss systemat- 
ics in the dust obscuration with stellar population age 
and bolometric and UV luminosity. F or ease o f com - 
parison with other studies, we assume a ISalpeteH (|1955ft 
initial mass function (IMF) and adopt a cosmology with 
H = 70kms- 1 Mpc" 1 , fl A = 0.7, and Q m = 0.3. 

2. SAMPLE SELECTION AND HERSCHEL PACS 100 AND 
160 /zM DATA 

2.1. Rest-Frame UV- Selected Sample of Star- Forming 
Galaxies at z ~ 2 — 3 

The GOODS-North field was targeted as part of an on- 
going imaging and spe ctroscopic survey o f UV-selected 
galaxies at z ~ 2 — 3 (jSteidel et all 120041 primarily to 
take advantage of the extensive multi-wavelength data 
that exist for this field. Spectroscopic catalogs and anal- 
ysis of the GOODS-North sample, including the stellar 
masses and dust attenuation of L* g alaxies at z ~ 2 — 3, 
are presented in IReddv et all ((2006a). Briefly, the galax- 
ies were selected to lie at redshifts 1.5 < z < 3.4 based 
on the "BM," "BX," and Lyman Break galaxy (LBG) 
color criteria (jSteidel et al.l f2003. 2001 lAdelberger et all 
2004). Extensive spectroscopic followup of candidates 
brighter than 1Z = 25.5 was conducted using the blue 
arm of the Low Re solution Imaging Spectrograph (LRIS; 
ISteidel et"aT] |2004') on the Keck I telescope. 

The criteria used to construct o ur sample are identi- 
cal to those of IReddv et al.l ()2010ft . Only galaxies with 
secure spectroscopic redshifts 1.5 < z spoc < 2.6 were con- 
sidered as this is the redshift range over which the MIPS 
band is sensitive to rest- frame 8/xm emission. Galax- 
ies with any AGN signatures based on optical emission 
lines (Lya, CIV, NV) or a power-law SED through the 
IRAC and MIPS 24 /an bands were excluded. Addition- 
ally, galaxies were included only if we were able to obtain 
robust point spread function (PSF) fits to their 24 /im 
emission; in practice, this meant that a few galaxies were 
rejected because they are co nfused with nearby neigh- 
bors (see Red dv et all l201fJI for the parameter used to 
measure the degeneracy of the PSF fits). This sample 
forms the basis of our stacking analysis, and includes a 
total of 146 galaxies, approximately 30% of which are de- 
tected with greater than 3 a significance in the GOODS- 
N 24 /xm data. 

2.2. Properties of the Sample and SED Modeling 

IReddy fe Steidell (|2009ft use a larger spectroscopic and 
photometric sample (the latter extending to 1Z ~ 27.0) 
in 31 independent fields to quantify the selection func- 
tion for the UV sample and determine the UV luminosity 



function at z ~ 2. This analysis indicates that UV- 
selected galaxies commonly targeted for spectroscopy 
with Tt < 25.5 have luminosities similar to L* of the 
UV luminosity functio n, where LJj V « 4 x 10 10 L Q 
(jReddv fc Steidelll2009f ). Throughout the text, L* refers 
to the characteristic luminosity of the UV luminosity 
function, unless stated otherwise. The stellar popula- 
tions of these galaxies have been modeled p reviously, and 
we refer the reader to iReddv et al.l (|2010l ) for a detailed 
description of the modeling procedure. 

Briefly, Chariot & Bruzual (in prep.) stellar popu- 
lation models with a range of exponentially-declining 
star formation histories r = 10, 20, 50, 100, 200, 500, 
1000, 2000, and 5000 Myr, as well as a constant star for- 
mation history, were fit to the observed U n G7Z, JK S , 
and Spitzer/IRAC data. In addition, we considered ages 
spaced roughly logarithmically between 70 and 5000 Myr, 
excluding any that are older than the age of the Universe 
at the redshift of each galaxy. The lower limit to the al- 
lowed ages (70 Myr) was adopted to reflect the typical 
dynamical time-scale as inferred from ve locity dispersion 
and size measur ements of z ~ 2 LBGs ()Erb et al.l 120061 : 
lLaw et aD 120071). Redden i ng wa s incorporated by em- 
ploying the ICalzetti et al.l (|2000l ) attenuation curve and 
allowing E(B — V) to range between 0.0 and 0.6. 

We adopt the best-fit stellar population parameters ob- 
tained when assuming a constant star formation (CSF) 
history, unless an exponentially declining (r) model gives 
a significantly better fit to the broadband data. Gener- 
ally, the x 2 values assuming the CSF model were sim- 
ilar to those obtained when r is allowed to vary freely. 
Further, more extreme star formation histories, where 
the ratio of the age to r is much greater than unity, 
i.e., t agc /r 3> 1, can be ruled out based on the presence 
of O star and Wolf-Rayet features in the composite UV 
spectra of z ~ 2 galaxies, as well as the fact that such 
models often predict ages that are unrealistically much 
younger than the dynamical timescale at z ~ 2. For the 
present analysis, we use the results of the stellar pop- 
ulation modeling to distinguish those galaxies that have 
young star formation ages < 100 Myr. Such galaxies have 
been shown previously to depart from the UV attenua- 
tion curve established for more typical galaxies at z ~ 2 
(jReddv et al.ll2010l l2006bl) . and we wish to explore this 
difference with the Herschel data. 

2.3. Construction of Subsamples 

The sample of 146 UV-selected galaxies was subdivided 
into subsamples in order to investigate the differences in 
dust attenuation between galaxies with blue and red UV 
spectral slopes, those with high bolometric luminosities 
(Lbol — ^m + £tjv); an d those with young stellar popu- 
lation ages. The properties of these subsamples, includ- 
ing the criteria used to construct them, the number of 
galaxies in each subsample, and their (3 and redshift dis- 
tributions, are summarized in Table [TJ As noted above, 
the ages are estimated from stellar population modeling 
of the rest-frame UV through Spitzer IRAC photometry. 
For the purpose of contructing the subsam ples, we esti- 
mated bolometric luminosities based on the lReddv et al.l 
(2010) calibration between 8/xm and infrared luminosity 
(see next section). Finally, UV slopes (3 were determined 
from the G — 1Z colors of galaxies as follows. We gen- 
erated power laws in /(A) ex A* 3 for —2.5 < (3 < 1.0 



with A/3 = 0.01. These w ere atte n uated for the Lya 
forest opacity assuming the iMadaul (|1995| ) prescription 
and multiplied by the G and 1Z transmission filters. The 
G band is affected by the Lya forest only for redshifts 
z = 2.5; statistical fluctuations in the forest will not af- 
fect (3 as most of the galaxies considered here lie below 
redshift z < 2.5. For this same reason, Lya lies outside 
of the G-band filter and will therefore not affect (3. The 
UV slope for a given galaxy is taken to be the one which 
gives the closest match in G — TZ color to the observed 
value. The error in UV slope is related directly to the 
error in color and is typically o~p ~ 0.11. 

2.4. Previous MIPS Results 

To provide a context for our present analysis, we 
briefly summarize previous efforts to constrain the dust 
emissio n and bolome t ric lum inosities of galaxies at z ~ 
2 — 3. IReddv et al.l (|2006bD investigated the variation 
in 24 fim. flux with dust-corrected UV and X-ray mea- 
sures of the bolometric luminosities of the same UV- 
selecte d galaxies analyzed here. Extending upon this 
result, IReddv et al.l (|2010l ) used a sample of 90 Lyman 
Break galaxies to examine the relationship between rest- 
frame 8/xm (z/L„[8/xm] = Lg) and Ha luminosity at 
z ~ 2, finding a tight trend between the two, with a 
scatter of k, 0.24 dex. Using Ha luminosity as a proxy 
for tot al star formation rat e - after accounting for dust 
with a ICalzetti et al.l (|2000h attenuation curve - then al- 
lowed these authors to establish a relationship between 
Lg, and dust obscured SFR, or Lir. 

The previous studies found that typical (£*) star- 
forming galaxies at z ~ 2 have dust attenuations, or 
infrared-to-UV luminosity ratios, Lir/Luv — 5, similar 
to the values predicted from the UV slope , /?, using the 
local correlation between (3 and -Lir/Luv (|Meurer et al.l 
1999). Investigations of the 24 (im emission of galaxies 
selected at rest-frame optical wavelengths - resulting in 
samples that are not orthogonal to the one analyzed here 
- have reached similar conclusions regarding the validity 
of UV-based dus t corrections for m oderately-luminous 
galaxies at z ~ 2 (|Daddi et al.ll2007fl . These resu l ts have 
been extended to higher redshift: iMagdis et al.l (|2010al) 
demonstrate that UV-based dust corrections for LBGs at 
z r~j 3 yield bolometric luminosities comparable to those 
inferred from infrared, radio, and millimeter measures. 
The 160 /jm emission from z ~ 3 LBGs also corresponds 
to Lib. that are similar to th ose obtained with 24 /im esti- 
mates (jMagdis et al.ll2010bh . While this is important val- 
idation of our understanding of star forma tion and dust 
obscuration at high redshift, the LBGs that Magdis et al. 
(2010b) detected in stacked 160 /im data have signifi- 
cantly higher luminosities ((£ir) ~ 1.6 x 10 12 L Q ) and 
star formation rates (> lOOM yr _1 ) than those of the 
typical, L* , galaxies analyzed here. 

In any case, while this agreement is encouraging, and 
instills confidence in our ability to recover the dust at- 
tenuation of typical galaxies from measurements of their 
UV continuum slopes, uncertainties in the fc-corrections 
and the conversion between 8 /im and infrared luminosity 
(as well as conversion between X-ray emission and SFR, 
for the X-ray stacking analyses) no doubt introduce some 
scatter in the inferred dust obscurations. The primary 
goal of our present analysis to obtain more direct mea- 
sures of the total infrared luminosities of z ~ 2 galaxies 



TABLE 1 
Properties of the Stacks I: 100 /xm Fluxes, Redshifts, and UV Slopes 



Sample 


Criterion a 


/loo 
rexc = 3'.'43 b 


/lOO 
No Exclusion c 




A. 

ALL 

UV- 

SELECTED 


All 


(3.0 ±0.5) x 10" 4 Jy 
N = 106 

(ft) = -1.41 ±0.42 
(z) = 2.08 ±0.24 


(3.0 ± 0.4) x 10" 4 Jy 
N = 146 

(8) = -1.41 ±0.41 
(z) =2.08 ±0.26 




B. 
L* 
GALAXIES 


Age > 100 Myr; 
L bol < 10 12 L 


(2.4 ±0.5) X 10~ 4 Jy 
N = 83 

IB) = -1.48 ±0.37 
(z) = 2.09 ±0.24 


(2.8 ±0.5) x 10~ 4 Jy 
N = 114 

(ft) = -1.46 ±0.38 
(z) = 2.09 ±0.25 




C. 

BLUE 

UV-SLOPES 


Age > 100 Myr; 
^bol < 10 Lq; 

8 < -1.4 


(1.9 ±0.7) x 10" 4 Jy 
JV = 51 

(8) = -1.61 ±0.21 
(z) = 2.05 ±0.24 


(2.1 ±0.7) x 10" 4 Jy 
N = 67 

(8) = -1.59 ±0.21 
(z) = 2.05 ±0.26 




D. 

RED 

UV-SLOPES 


Age > 100 Myr; 
L bol <10 12 L Q ; 
> -1-4 


(3.8 ±0.9) x 10" 4 Jy 
N = 32 

(8) = -1.10 ±0.26 
(z) = 2.19 ±0.22 


(3.4 ±0.6) x 10" 4 Jy 
TV = 47 

(8) = -1.15 ±0.26 
(2) = 2.14 ±0.24 




E. 
ULIRGs 


Age > 100 Myr; 
L bol > 10 12 L 


(11.0 ±1.7) x 10" 4 Jy 

JV = 9 

{8) = -0.84 ±0.48 

<z} = 2.03 ±0.27 

<L bol } = (1.3±0.1) xlO 12 L 

L™ = 1.8 x 10 12 L 


(9.9 ±1.2) X 10~ 4 Jy 

N = 12 

(ft) = -0.84 ±0.46 

(2) = 2.22 ±0.29 

(i bo) ) = (1.3 ± 0.1) X10- 

L™ = 1.8 x 1O 12 L 


-4 T 

Lq 


F. 

YOUNG 

GALAXIES 


Age < 100 Myr 


3a: < 4.4 X 10~ 4 Jy 
JV = 14 

(8) = -1.17 ±0.40 
(z) = 2.00 ±0.26 


3ct: < 3.2 x HT 4 Jy 
N = 20 

(ft) = -1.30 ±0.40 
(2) = 2.00 ±0.28 





NOTE. — Each entry includes (a) the stacked flux at 100 /^m and its measurement uncertainty; (b) the 
number of galaxies contributing to the stack: and (c) the mean and sample dispersion of the UV slopes (/3) 
and redshifts of those galaxies. We assign 3 a upper limits to fluxes in cases where the 1 a measurement 
uncertainty is larger than the stacked flux. For Sample E (ULIRGs), we include the mean and error in 
the mean of the bolometric luminosity (Lbol! derived using the calibration of Redely ct al. 2010) and the 
maximum bolometric luminosity (L™^ x ) contributing to the stack. 

a The criteria used to select each subsamplc arc listed here. The ages were determined from SED-fitting 
to the broadband photometry (see text) and the bolometric luminosities, Lbol, were determined from our 
previous calibration between 24 (im and infrared luminosity dReddv et al.|[2010D . 

Criteria for excluding galaxies from the stack. If any galaxy lies within a distance r cxc of a nearby optical, 
i^ s -band, or IRAC source, it is excluded from the stack. 
c No exclusion radius adopted; all galaxies arc median stacked. 



using Herschel data, which are described below. 

2.5. Herschel and VLA Data 

The GOODS- Herschel Open Time Key Program (PI: 
D. Elbaz) includes « 124 hrs of 100 and 160 //m PACS 
imaging and 31 hrs of SPIRE imaging in the GOODS- 
North field. The 3 a depths of the PACS 100 and 160 /im 
images are 1.1 mJy and 2.6 mJy, respectively. None of 
the galaxies in the sample analyzed here are detected 
to these de pths. Furthe r deta ils on the data reduction 
are given in lElbaz et al.l f|2011h . PSFs were constructed 
based on the catal ogs of Herschel detections discussed in 
lElbaz et al.l ([201 1[ ). When performing PSF photometry, 
we adjusted all fluxes upward by multiplicative factors of 
1.37 and 1.29 at 100 and 160 /im respectively, to account 
for missing flux in the wings of the PSFs. The similar 
factor for the MIPS 24 ^m PSF is 1.22. 

The VLA radio 1.4 GHz data are described in 
iMorrison et al.1 (|2010T ) . Brief ly, a total of 165 hrs of obser- 
vations, including 42 hr from lRichard s (2000), were com- 
bined to produce a radio map of the GOODS-North field. 



The noise level of the image is ~ 3.9 /iJy beam -1 near the 
center and ~ 8/iJybeam -1 at 15' from the center. The 
synthesized beamsize of the radio data is ~ 1'.'7 x 1'.'6, 
corresponding to a physical scale of 14.1 x 13.3 kpc at 
z = 2.30. The radio map used in the stacking analysis is 
primary beam corrected. 

3. STACKING METHOD AND SIMULATIONS 

3.1. General Stacking Procedure 

The method used to stack the Herschel 100 and 
160 /im, Spitzer 24 /im, and the VLA radio data is as fol- 
lows. We extracted an area 50" x 50" around each object 
to be stacked. The area was chosen to be large enough to 
obtain a reliable estimate of the local background. These 
cutouts were then median combined after rotating each 
sub-image by 90° relative the previous sub-image con- 
tributing to the stack, in order to minimize image defects 
that are aligned with the scanning orientation of the data 
acquisition. 17 We adopted the median stacks to ensure 

17 To test for any systematic effect associated with asymmetries 
in the PSF, we also stacked without rotating each sub-image. The 



F(24|im] = 
27 ± 2 pJy 




F(lDO(irrl}= 
280±50(iJy 



F(160nm}= , 
590 ± 150 [Jjy 



S(1.4GHi)= 

2.0±0.4MJy 




Fig. 1. — PSF-convolved median stacked images for "Sample B" 
(L* galaxies) at 24, 100, 160 (im, and 1.4 GHz, with pixel scales of 
\'.'2, l'/2, 2'.'4, and 0?5, respectively. 

the results are not biased by bright outliers. We stacked 
the 24 /im data in exactly the same way as was done 
for the 100 and 160 /im data, irrespective of whether a 
galaxy was directly detected at 24 /im, in order to ensure 
the most consistent results across wavelengths. 18 

To test the effects of a non-negligible amount of flux 
in the wings of the stacked signal relative to the PSF 
(due to positional uncertainty), we performed both PSF- 
fitting photometry and simple aperture photometry on 
the stacked flux. The background and measurement un- 
certainty in the stacked image are taken respectively to 
be the mean and 1 a dispersion in flux computed by fit- 
ting many PSFs at random positions in the stacked im- 
age. We simulatenously fit the target flux (with a PSF) 
and a constant value to account for the background level 
in the stacked image. We found no significant system- 
atic bias in the aperture versus PSF-derived flux relative 
to the measurement uncertainties, and we adopted the 
PSF-determined flux measurements. For the radio data, 
the PSF-measured flux and total integrated flux differ by 
~ 20 — 30% due to bandwidth smearing. To account for 
this effect, we have adopted the fluxes computed from 
fitting the best-fit elliptical Gaussian profile (with peak 
intensity normalized to unity) to the stacked radio emis- 
sion, yielding the integrated flux density. 

3.2. Residual Images 

Residual maps at 100 and 160 /im were constructed in 
order to ensure the most robust determinations of the 
stacked fluxes. 19 Sources detected at > 3 a were sub- 
tracted from the science mosaics (none of the galaxies in 
our sample are detected at 100 and 160 /im at this level). 
We considered detections based both on a blind catalog 

results obtained with and without rotating the sub-images were 
identical within the uncertainties of the stacked flux measurements. 

18 For the most conservative estimate of the sta cked radio flux, 
we did not correct for bandwidth smearing (e.g., IPannella et al.l 
2009) , rather we fit an elliptical Gaussian to the stacked light profile 
and adopted the peak value for the radio flux. 

19 We did not construct a residual map for the VLA data given 
the higher resolution of these data and the lower 1.4 GHz source 
surface density. 



("blind-subtracted"), where a detection algorithm was 
used directly on the PACS mosaics, as well as a cata- 
log constructed using 24 /im priors to define the posi- 
tions of sources in the longer wavelength PACS bands 
( "prior-substracted" ) . In both cases, we used PSF pho- 
tometry to determine the fluxes of sources and subtract 
them from the science images. To ensure that there are 
no systematic effects that may bias the flux measure- 
ment made on the residual maps, we performed the same 
stacking analysis on the science images themselves. Com- 
parison between the stacks on the blind-residual maps, 
prior- residual maps, and the science images, showed that 
the stacked fluxes were within 10% of each other. Typi- 
cally, the stacked fluxes measured on the science images 
were lower than those measured on the residual images. 
This effect is attributed to the higher background level in 
the stack derived from the science images due to the PSF 
wings of nearby sources. The similarity in flux regardless 
of the image used for the stack suggests that the median 
combination of sub-images is robust to bright outliers, 
as expected. The flux measurements obtained by stack- 
ing on the prior-fit residual maps are adopted for further 
analysis. The median stacked data for L* galaxies (Sam- 
ple B; Tabled]) at 24, 100, and 160 /im, and 1.4 GHz, are 
shown in Figure [T] 

3.3. Exclusion Tests and Simulations 

A common concern in stacking data with relatively 
poor resolution is the potential contribution from un- 
related sources within the beam that are formally unde- 
tected (i.e., with a < 3 a significance) and which there- 
fore remain in the residual images, but which may con- 
tribute significantly to the stacked flux. If these unre- 
lated sources have a random spatial distribution, then 
their contribution to the stacked flux will be accounted 
for when we subtract the background in the stacked im- 
ages. However, if the unrelated sources are clustered with 
respect to the UV-selected galaxies, they will contribute 
to the stacked fluxes. To test for the effect of clustering 
on the stacked fluxes, we considered a number of tests, 
as we describe below. 

3.3.1. Exclusion Radius 

The most conservative measure of the median flux can 
be achieved by stacking only those galaxies without any 
nearby sources. The full-width-half-max (FWHM) of the 
PSF at 100 /im is small enough (~ 6") such that we can 
adopt an "exclusion radius" r oxc =FWHM/2, and still 
have enough galaxies to stack. We first constructed a 
catalog containing all optical, -Kg-band, and IRAC de- 
tections in the GOODS-North field. The 3er sensitivi- 
ties, as measured in a 2" diameter aperture in the opti- 
cal and Kg-band images, are 1Z ~ 27.6 and K s ~ 24.55. 
The 3 a sensitivities of the GOODS-North IRAC data, 
as measured in a 4" diameter aperture, are 26.56, 26.00, 
24.17, and 24.16, for the 4 IRAC channels, respectively. 
Any galaxies in our sample that lie within a distance 
r GXC of any source detected in the optical, near-IR, or 
with IRAC, are excluded from the stack. The sample 
itself was selected such that galaxies that were confused 
with any nearby MIPS 24 /im sources were excluded. 

The median 100 /im fluxes obtained with and with- 
out adopting an exclusion radius are listed in Table [TJ 
The differences between these fluxes are smaller than 



TABLE 2 
Properties of the Stacks II: Summary of Mid-IR, IR, and Radio Fluxes 





fli 


/ioo 


/l60 


/l.4 


Sample 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


A. 


(33 ±4) x 1CT 6 


(3.0 ±0.4) x 10~ 4 


(7.3 ±1.2) x 10~ 4 


(2.1 ±0.5) x 10 -6 


B. 


(27 ±2) x 10~ 6 


(2.8 ±0.5) x 10~ 4 


(5.9 ±1.5) x 10~ 4 


(2.0 ±0.4) x 10~ 6 


C. 


(21 ±3) x 10" 6 


(2.1 ±0.7) x 10~ 4 


(5.4 ±1.6) x 10~ 4 


(1.6 ±0.6) x 10~ 6 


D. 


(35 ±3) x 10" 6 


(3.4 ±0.6) x 10~ 4 


(6.6 ±2.1) x 10~ 4 


(2.6 ±0.7) x 10~ 6 


E. 


(130 ±4) x 10" 6 


(9.9 ± 1.2) x 10~ 4 


(22.4 ±3.8) x 10~ 4 


(6.6 ±1.3) x 10~ 6 


F. 


(10 ±2) x 10 -6 


3<r: < 3.2 X 10~ 4 


(4.6 ±3.3) x 10~ 4 


3ct: < 3.1 X 10~ 6 



NOTE. — The quoted errors reflect measurement uncertainty in the stacked fluxes. 



TABLE 3 
Properties of the Stacks III: Biases and Errors in Stacked Fluxes 





24 fira 


100 (im 


160 /an 


1.4 GHz 


Sample 


1-Bias a / Error b 


1-Bias a / Error b 


1-Bias a / Error b 


1-Bias a / Error b 


A. 


0.94 / 0.002 


0.96 / 0.012 


0.92 / 0.014 


1.03 / 0.016 


B. 


0.94 / 0.003 


0.95 / 0.017 


0.92 / 0.022 


1.02 / 0.024 


C. 


0.94 / 0.015 


0.96 / 0.039 


0.92 / 0.040 


1.03 / 0.054 


D. 


0.94 / 0.004 


0.96 / 0.035 


0.92 / 0.045 


1.03 / 0.038 


E. 


0.95 / 0.009 


0.95 / 0.046 


0.93 / 0.058 


1.05 / 0.075 



a Average bias of stacked flux, defined as the ratio of the mean measured flux to the simulated flux, 

(r" s )/f m . 

Fractional uncertainty in the mean stacked flux, taken as the ratio of the error in the mean of the 
measured stacked flux (er/\/]V) to the median measured flux, i.e., cr(/ 771eas )/(v / ^V(/ meas )). 



the 1 a measurement uncertainties in the stacked fluxes. 
The similarity in flux between the exclusion and no- 
exclusion cases implies that any objects that may clus- 
ter around the UV-selected galaxies do not contribute 
significantly to the stacked far-infrared fluxes of the UV- 
selected galaxies. We note that there may exist very faint 
sources that are undetected in the optical, near-IR, and 
IRAC data and which may lie close to our targets. How- 
ever, we consider it unlikely for sources that are faint 
at virtually all other wavelengths (optical, near-IR, and 
the Spitzer IRAC and MIPS bands) to be bright enough 
at 100 /mi to contribute significantly to the stacked flux 
of the UV-selected galaxies. Performing a similar test 
at 160 /im is not possible, as the larger FWHM results 
in an exclusion radius that precludes any galaxies from 
being stacked. Because the 100 and 160 pirn emission 
arises from the same mechanism (dust emission), then 
one would conclude that the effects of clustering are un- 
likely to be significant at 160 /im if the 100 pirn stacked 
emission yields similar infrared luminosities to those ob- 
tained from the stacked 160 /im data. In Section |4] we 
present evidence that suggests that an additional con- 
tribution to the 160 /im flux from other sources must 
be negligible compared to the flux from the UV-selected 
galaxies of interest here. Based on this evidence, and to 
take advantage of all the galaxies in our sample, we pro- 
ceed by adopting the fluxes obtained without using the 
exclusion radius; these values are given in Table [2] Be- 
low, we discuss two additional tests used to verify these 
stacked fluxes. 

3.3.2. Comparison with Random Stacks 

The probability of a chance measurement that results 
in fluxes as high as the ones obtained by stacking on the 
positions of UV-selected galaxies (i.e., the target stacks) 



can be determined by stacking at random positions in 
the images (i.e., random stacks). We performed this ran- 
dom stack test on the radio and 24, 100, and 160 /mi data 
(we did not exclude positions that correspond to detected 
sources). The results for the mid-IR and IR simulations 
for Sample B (stacking on N = 114 positions) are shown 
in Figure [2] For this sample, the 24 and 100 /im fluxes 
measured for 10,000 random stacks were never as high 
as those obtained for the target stacks. For the 160 /mi 
data, 15 out of 10,000 random stacks result in fluxes that 
were within the 1 a measurement uncertainty of the tar- 
get stack. The random stacks indicate a very low proba- 
bility (< 0.2%) of accidentally recovering stacked fluxes 
as high as the ones observed for the sample of L* galax- 
ies at z ~ 2, and these probabilities are consistent with 
those expected based on the measurement uncertainties 
of the stacked fluxes. 

3.3.3. Stacks of Simulated Galaxies 

Building on the random stack tests, point sources of 
known flux density were added at random locations in the 
residual images (using the same PSFs that are used to 
obtain photometry 20 ) and were recovered by a stacking 
analysis. We added the same number of point sources to 
the residual images as are used in constructing the stacks 

20 We performed another test to determine if asymmetries in the 
intrinsic Herschel PSF result in systematic differences in the pho- 
tometry obtained for simulated sources. For this test, we redid the 
simulations by adding point sources of known flux density, where 
we assumed the flux profile given by the PACS PSF measured by 
observing the asteroid Vesta. The sources were then recovered us- 
ing the PSF measured from the GOODS-N Herschel images. Based 
on this second set of simulations, we find less than 5% systematic 
offset between the difference of input and output flux relative to the 
difference obtained when assuming the same PSF, both for adding 
sources to the images and recovering their fluxes. 
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Fig. 2. — Comparison of measured stacked fluxes (mean and uncertainty indicated by the solid lines and hashed regions, respectively) 
for Sample B and the distribution of fluxes obtained by stacking on N = 114 random positions, repeated 10000 times (histograms). 
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Fig. 3. — Distribution of stacked fluxes for artificial galaxies added to the 24, 100, and 160 (im images (histograms) for Samples A through 
E, relative to the simulated flux (numbers in panels). 



for Samples A through E. Flux densities were assigned 
based on a Gaussian distribution with a mean equivalent 
to the median stacked flux for each Sample (the stacks 
are insensitive to the dispersion in simulated fluxes) . The 
ratio of the measured stack to simulated mean flux are 
shown in Figure [3] for the different samples, with num- 
bers and input mean fluxes indicated, for the 24, 100, and 
160 /xm data. This ratio is close to, but not exactly, equal 
to 1, with a bias of just a few percent. The biases and 
dispersions measured for the different stacks are summa- 
rized in Table [3] The error in the mean recovered stacked 
flux from the simulations is much smaller than the mea- 
surement uncertainty. Therefore, for the subsequent dis- 
cussion we assume a flux uncertainty that is equal to the 
measurement uncertainty of the stacked flux. Finally, we 
corrected the observed fluxes (listed in Table ^ by the 
bias factors given in Table [3J before inferring the total 
infrared luminosities, as we discuss in the next section. 



3.3.4. Summary of Tests 

We have performed several tests to validate our mea- 
sures of the stacked fluxes for UV-selected galaxies. The 
method employed here is able to recover the stacked 
fluxes of galaxies with a bias of just a few percent rela- 
tive to the known fluxes of artificial sources added to the 
1.4 GHz, and 24, 100, and 160 /xm images. These tests 
also imply a very low probability (< 0.2%) of accidentally 
recovering stacked fluxes that are as high as the target 
stacked fluxes. Finally, we test for the effect of source 
clustering on the stacked flux at 100 /im. By excluding 
galaxies from the stack which have any nearby sources, 
we find that the remaining galaxies have a median flux 
that is identical within the measurement errors of the 
median flux inferred for the entire sample of galaxies. 
This indicates that any sources that may cluster around 
the UV-selected galaxies (if they exist) do not affect the 
stacking results. 



4. INFRARED LUMINOSITIES AND DUST SEDS 

4.1. Dust SED Templates and Extrapolation to the 

Radio 

To infer the total infrared luminosities L(8 — 
1000 /im) = Lir of the z ~ 2 galaxies, we employed 
several publicly availabl e dust SED templ ates including 
those of ICharv fc Elbad (12001 C'C EOl" ).IDale fc Heloul 
(|200l f"DH02"). and lRieke et all (|200l ("Rieke+09"). 
We also in c lude the median templates presented in 
lElbaz et al.l (|2011l ). Specifically, these authors define an 
"infrared main sequence" of galaxies, where the ratio of 
Lir to Lg (IR8) is universal for most star-forming galax- 
ies at z < 2.0, having a value of Lir/Ls = 4.9^ 2 |. The 
"infrared starbursts" are considered to be those galax- 
ies with IR8 ratios in excess of « 15. El baz et al.l (2011) 
demonstrate that variations in IR8 may be due primarily 
to differences in the infrared luminosity surface density 
of galaxies, such that main sequence galaxies are char- 
acterized by star formation that is more extended than 
that present in starbursts (see Section [5]). We fit the me- 
dian IR SED of these main sequence ("Elbaz+11-MS") 
and starburst galaxies ("Elbaz+11-SB") to the stacked 
fluxes. 

The lReddv et al.l ()2010f ) calibration between rest-frame 
8/im emission and infrared luminosity is considered as 
well. In this case, Lg, is computed by k-correcting the 
24 /im flux using the average of 12 loc al galaxy mid-IR 
SEDs specified in lReddy etaU (|2006bh . The correlation 
between L% and Ha luminosity is then used to infer the 
Lg-Lm conversion, based on using Ha as an indepen - 
dent probe of the star formation (|Reddv et al.l [2010). 
Finally, we make use of the radio-infrared correlation to 
provide an independent estimate of Lir, assuming that 
this relation does not evolve with redshift, as is consis- 
tent with the evidenc e , at least at redshi f ts z < 3 (e.g., 
Appleton et al.l 12004 llvison et al.l I2010J: ISargent et alj 
20101: iBournc et al. 20111: iMao et al.ll2011ft. For the sub- 



sequent ana lysis , we u se the radio-in frared correlation 
published in[Bei3 (|2003ft . The sample of lYun et al.l (|200lD 
is roughly a f actor of 10 larger than the one analyzed 
bvlBell (|2003f ). However. lYun et al.l (|2001l ) calibrate the 
60 /im luminosity, Lqq, with radio luminosity, and the 
former requires some assumption about the relation be- 
tween Lqq and Lir. To make the min imum numb er of 
assumptions, we therefore adopted the Bell (2003) cali- 
bration which directly relates the total infrared luminosi- 
ties to the specific luminosity at 1.4 GHz. 

The CE01 and Rieke+09 models are parameterized 
such that we can relate any template to a given total in- 
frared luminosity. The Elbaz+11-MS/SB templates are 
normalized to an infrared luminosity of Lir = 10 n L0. 
The DH02 models are presented as a function of radi- 
ation field intensity and, following the literature (e.g., 
iPapovich et al.l 120071 ). we associate infrared luminosities 
with each of the DH02 template s assu ming the empiri- 
cal calibration of iMarcillac et al.1 (|2006ft . The CE01 and 
Rieke+09 models also include an extrapolation to ra- 
dio wavelengths by assuming some version of the local 
radio - far-infrare d correlation ([Condon! 119921 : lYun et all 
12001b iBelll I2003T ) . For a consistent treatment, we have 
adjusted the models (or added to them, in the case of 
DH02 and Elbaz +11-MS/SB) fo r a radio spectrum with 
index 7 = —0.8 (|Condonl ll992). normalized to have a 



specific luminosity at 1.4 GHz c orrespondin g to the Lir 
for that template, assuming the Bell (2003) calibration. 

4.2. Fitting Procedure 

Because the far-infrared peak of the dust emission 
in the SED templates shifts to shorter wavelengths at 
higher luminosities, we can either (1) fit these tem- 
plates directly to the observed infrared and radio fluxes 
( "luminosity- matched" fitting), or (2) find the template 
that best matches the infrared colors and then normal- 
ize this template to the observed fluxes ( "color-matched" 
fitting). The color-matched fitting yields larger uncer- 
tainties in Lir because the color errors are larger than 
those for individual flux measurements. We have adopted 
both methods for comparison purposes. The Elbaz+11- 
MS/SB templates are based on the composite IR spectral 
energy distribution of the "main sequ e nce" and "star- 
burst" galaxy samples of lElbaz et al.l (|2011h . and are 
normalized to Lir = 10 11 L Q . We assume the same spec- 
tral shape (i.e., no color dependence) when fitting these 
templates to the observed fluxes. Further, in finding the 
best-fit template, we include results where all fluxes have 
been weighted equally, as well as being weighted by their 
errors. The equal weighting is done to ensure that the 
higher S/N 24 /im data (and thus the smaller measure- 
ment uncertainties at 24 /im) do not unduly skew the 
template fits. 

4.3. Infrared Luminosities 

Table |4] lists the I/ir corresponding to the best-fit tem- 
plates for different combinations of the stacked fluxes 
for Samples A through F, after taking into account the 
biases and dispersions in these fluxes (see discussion 
above and Table [3]) . Fitting for different combinations 
of the infrared and radio data allows us to determine 
if the inclusion of any of the stacked fluxes results in 
significant changes in the inferred infrared luminosities. 
Als o included in this table are the results from using 
the lReddv et al.l (|2010t ) calibration between L 8 and Lir, 
and t he Li r corresponding radio luminosity assuming the 
iBelll ([2003) radio-IR correlation. 

The systematic uncertainties in log(LiR) for the differ- 
ent combinations of templates, fluxes, and fitting meth- 
ods for each of the subsamples are typically within 
0.1 dex. We note that the Lir derived from 160 /im is 
not systematically larger than that inferred from 100 /im, 
and hence runs counter to what one might expect if 
the larger 160 /zm beam included more sources clustered 
around the UV-selected galaxies that contribute to the 
measured far-infrared flux. The stacked radio data have 
significantly higher resolution (~ l'/7 at 1.4 GHz versus 
~ 11" at 160 /im) and yield a radio flux that implies an 
Lir that is also not any smaller than the value inferred 
from the 160 /im data. Based on this evidence, we con- 
clude that the 160 /im emission from sources proximate 
to the UV-selected galaxies is negligible compared to the 
emission from the UV-selected galaxies themselves. 

The notable outliers listed in Table 0] are the Lir de- 
termined from 24 /im flux density alone, which tend to 
overpredict Lir, particularly for the highest luminosity 
subsamples D and E, relative to cases where we include 
the 100 and 160 /im data in most of the template fits 
(Figure 0}. Alternatively, with the Elbaz+11-MS tem- 
plate, 24 /zm-only determinations result in systematically 



TABLE 4 
Comparison of Infrared Luminosities (£ir) 



Template 



A, v 



Sample A Sample B Sample C Sample D Sample E Sample F 



R10a c 
Bell03 d 

Elbaz+11-MS Lum c 



Elbaz+11-MS Lum-Eq Weight f 
Elbaz+11-SB Lum e 



Elbaz+11-SB Lum-Eq Weight f 
CE01 Lum e 



CE01 Lum-Eq Weight f 
CE01 Col s 

CE01 Col-Eq Weight h 

DH02 Lum e 



DH02 Lum-Eq Weight f 
DH02 Col s 

DH02 Col-Eq Weight h 

Rieke+09 Lum c 



Rieke+09 Lum-Eq Weight f 
Rieke+09 Col s 

Rieke+09 Col-Eq Weight h 



24 

1.1 

24 

100 

160 

100, 160 

24, 100, 160 

24, 100, 160, 1.4 

24, 100, 160, 1.4 

24 

100 

160 

100, 160 

24, 100, 160 

24, 100, 160, 1.4 

24, 100, 160, 1.4 

24 

100 
160 

100, 160 
24, 100, 160 
24, 100, 160, 1.4 
24, 100, 160, 1.4 
100, 160 
24, 100, 160 
24, 100, 160, 1.4 

24 

100 
160 

100, 160 
24, 100, 160 
24, 100, 160, 1.4 
24, 100, 160, 1.4 
100, 160 
24, 100, 160 
24, 100, 160, 1.4 

24 

100 
160 

100, 160 
24, 100, 160 
24, 100, 160, 1.4 
24, 100, 160, 1.4 
100, 160 
24, 100, 160 
24, 100, 160, 1.4 



2.0 + 0.3 

2.2 + 0.5 

1.9 + 0.2 
3.0 + 0.4 
3.0 + 0.4 
3.0 + 0.3 

2.3 + 0.3 

2.3 + 0.3 

3.0 + 0.2 

3.7 + 0.5 

2.4 + 0.3 

2.3 + 0.3 

2.4 + 0.2 
2.6 + 0.3 

2.6 + 0.4 

2.3 + 0.2 

3.4 + 0.4 
2.4 + 0.3 

2.2 + 0.3 

2.3 + 0.3 

2.4 + 0.3 
2.4 + 0.4 

2.2 + 0.2 

2.1 + 0.4 

2.3 + 0.6 

2.4 + 0.7 

3.0 + 0.4 

2.1 + 0.3 

2.2 + 0.3 

2.1 + 0.2 

2.3 + 0.3 

2.3 + 0.3 

2.2 + 0.2 

2.2 + 0.4 

2.4 + 0.7 

2.3 + 0.6 

5.0 + 0.6 

2.5 + 0.3 

2.2 + 0.3 

2.3 + 0.2 
2.3 + 0.2 
2.3 + 0.3 
2.2 + 0.2 

3.1 + 0.6 

2.7 + 0.7 

2.6 + 0.7 



1.8 + 0.1 

2.2 + 0.4 

1.6 + 0.1 

2.9 + 0.5 

2.4 + 0.6 

2.7 + 0.4 
1.7 + 0.2 
1.7 + 0.2 

2.5 + 0.2 

3.1 + 0.2 

2.3 + 0.4 
1.9 + 0.5 

2.1 + 0.3 

2.7 + 0.3 

2.6 + 0.3 
2.0 + 0.1 

2.6 + 0.2 
2.3 + 0.4 

1.8 + 0.4 
2.0 + 0.3 
2.3 + 0.3 
2.3 + 0.3 
1.8 + 0.1 
1.8 + 0.5 

2.2 + 0.8 

2.2 + 0.8 

2.3 + 0.2 

2.0 + 0.4 

1.8 + 0.4 

1.9 + 0.2 
2.2 + 0.2 

2.2 + 0.3 

1.8 + 0.1 

1.5 + 0.5 

2.1 + 0.7 

2.2 + 0.8 

3.6 + 0.3 

2.3 + 0.4 

1.9 + 0.4 
2.1 + 0.2 

2.3 + 0.3 

2.3 + 0.3 
1.9 + 0.1 
2.8 + 0.8 

2.4 + 0.8 
2.3 + 0.8 



1.2 + 0.2 
1.7 + 0.6 

1.2 + 0.1 

2.0 + 0.6 

2.1 + 0.6 
2.1 + 0.4 

1.3 + 0.2 
1.3 + 0.2 
2.1 + 0.2 

2.3 + 0.3 

1.6 + 0.5 

1.7 + 0.5 

1.6 + 0.3 
2.0 + 0.3 

2.0 + 0.4 

1.7 + 0.2 

1.7 + 0.2 
1.6 + 0.5 

1.5 + 0.4 

1.6 + 0.3 
1.6 + 0.3 
1.6 + 0.4 

1.5 + 0.2 

2.1 + 1.0 

1.6 + 0.9 

1.6 + 0.9 

1.7 + 0.2 

1.5 + 0.5 

1.6 + 0.4 

1.5 + 0.3 

1.6 + 0.2 
1.6 + 0.3 

1.6 + 0.2 

1.7 + 0.7 

1.5 + 0.8 

1.6 + 0.9 

2.0 + 0.2 

1.8 + 0.6 

1.7 + 0.5 

1.7 + 0.2 

1.8 + 0.3 
1.8 + 0.4 

1.7 + 0.2 
2.4+1.3 

1.8 + 1.0 
1.8 + 0.9 



2.7 + 0.2 
2.9 + 0.8 

2.2 + 0.2 

3.8 + 0.7 

2.9 + 0.9 

3.5 + 0.5 

2.3 + 0.3 
2.3 + 0.3 
3.0 + 0.2 

4.3 + 0.4 

3.0 + 0.5 

2.3 + 0.7 
2.7 + 0.4 

3.6 + 0.4 

3.5 + 0.5 

2.4 + 0.2 

4.1 + 0.3 
3.1 + 0.6 

2.1 + 0.6 

2.7 + 0.4 
3.3 + 0.4 

3.2 + 0.4 

2.3 + 0.2 
1.7 + 0.6 
3.0 + 1.2 

3.0 + 1.2 

3.6 + 0.3 

2.7 + 0.5 

2.1 + 0.6 

2.5 + 0.4 
3.0 + 0.4 

3.0 + 0.5 

2.2 + 0.2 

1.7 + 0.6 

2.8 + 1.1 

3.6 + 1.5 

6.1 + 0.5 

3.0 + 0.5 

2.1 + 0.6 

2.5 + 0.3 

3.1 + 0.4 

3.0 + 0.5 

2.2 + 0.2 

3.6 + 1.3 
3.2 + 1.3 

3.1 + 1.2 



15.7 + 0.5 0.3 + 0.1 



7.9 + 1.7 3a: < 3.0 



9.0 + 0.3 
12.7 + 1.4 
11.1 + 1.7 
12.0+ 1.0 

9.2 + 0.4 

9.2 + 0.5 

11.3 + 0.3 

17.7 + 0.5 

10.0+ 1.1 

8.6+ 1.3 

9.4 + 0.9 

15.5 + 0.7 

15.0 + 0.8 
8.8 + 0.2 

41.3+ 1.2 

10.9+ 1.2 

7.9+ 1.2 

9.4 + 0.9 

12.3 + 0.4 

12.1 + 0.6 
8.2 + 0.2 
7.9+ 1.5 

10.4 + 2.3 

12.1 + 2.7 

21.5 + 0.6 
7.8 + 0.9 
7.4+ 1.2 
7.7 + 0.6 

13.6 + 0.7 

13.2 + 0.8 

7.4 + 0.2 
6.9+ 1.3 

10.5 + 2.4 

9.5 + 2.2 

46.1 + 1.3 

8.2 + 0.9 

6.1 + 0.9 
7.1 + 0.6 

11.2 + 0.7 

11.0 + 0.9 

6.3 + 0.2 
12.5 + 2.7 
11.9 + 2.7 

12.1 + 2.7 



0.5 + 0.1 
3<r: < 2.8 

1.7 + 1.1 
3cr: < 2.8 
3<r: < 2.8 
3<r: < 2.8 
3cr: < 2.8 

1.0 + 0.2 
3cr: < 2.2 

1.3 + 0.9 
3<r: < 2.2 
3cr: < 2.2 
3<r: < 2.2 
3cr: < 2.2 

0.4 + 0.1 
3<r: < 2.3 

1.2 + 0.8 
3cr: < 2.3 
3<r: < 2.3 
3<r: < 2.3 
3cr: < 2.3 



0.6 + 0.1 
3<r: < 2.0 

1.3 + 0.9 
3cr: < 2.0 
3<r: < 2.0 
3<r: < 2.0 
3cr: < 2.0 



0.5 + 0.1 
3<r: < 2.3 

1.4 + 1.0 
3cr: < 2.3 
3<r: < 2.3 
3cr: < 2.3 
3cr: < 2.3 



a Luminosities arc in units of 10 Lq. Errors in luminosities arc derived from the uncertainties in the stacked fluxes (see text). 

Wavelengths (in micron) or frequency (in GHz) used to compute the best-fit Lir. 
c Lir determined from the calibration of 24 //m with Lir of Roddy ct al. (2010), which is based upon the correlation between 8 [im and Ha 
luminosity for L* galaxies at z ~ 2. 

Lir determined from the radio-IR correlation of Bell (2003), and assuming a radio spectral slope of 7 — —0.8. 
c Lir for the lElbaz et all J2011) main sequence /starburst. ICharv &: Elbazl I200U) , lDale &: HeloTJ 120021) , or IRieke et al.l 1 12009) template that 
best matches the observed fluxes (weighted by the flux errors), irrespective of the infrared colors. 

f Lir for the lElbaz et al.l 120111) main sequence/starburst. lCharv fe Elbazl (f200ll) . IDale fc Heloul | [2002D . or lRieke et al.l 120091) template that best 
matches the observed fluxes (including 1.4GHz), with all fluxes given equal weight in the fitting, irrespective of the infrared colors. 
s Lir for the Chary & Elbaz (2001), Dale ik Hclou (2002), or Rickc ct al. (2009) template that best fits the observed infrared colors and is 
normalized to match the observed fluxes, with colors and fluxes weighted by their errors. 

11 Lir for the Chary & Elbaz (2001), Dale & Hclou (2002), or Rickc ct al. (2009) template that best fits the observed infrared colors and is 
normalized to match the observed fluxes (including 1.4 GHz), with colors and fluxes weighted equally. 
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Fig. 4. — Comparison between the Lir computed from the 24, 
100, 160 (im, and 1.4 GHz equally-weighted fluxes, and the Ljr 
computed from the 24 (im d ata only. The h eavy black points denote 
the Ljr computed from the Rcddy et al. (2010) calibration versus 
the Ltr, obtained from the best-fit DH02 template that includes all 
the data (24, 100, 160 (im, and 1.4GHz), equally-weighted. Error 
bars reflect the uncertainty in the Ljr inferred from the stacked 
measurements, except for the high luminosity subsample (Sample 
E), where the errors include the intrinsic dispersion in 24 /im fluxes 
of objects that contribute to the sample. For clarity, error bars are 
not shown for the Elbaz+11-MS/SB points. 

lower Lir. These differences are due to the intrinsic vari- 
ation in the ratio of Lg to Lir luminosity present in 
the templates relative to the observed SED. The system- 
atic overestimation of Ltr based on Lg alone has been 
noted before fo r ultraluminous infrared galaxies stud- 
ied with Spitzer (|Papovich et al.ll2007t iDaddi e t al. 2007; 
Papqvich et al. 2009t iMurphv et al.ll2011b iM agnelli et al. 
20111) and Herschel i^ovdon et al.l 120101: lElbaz et al.l 
20101) . For the samples considered here, the same 24 /im 
bias is apparent for all templates (except for Elbaz+11- 
MS) even at LIRG luminosities, where the difference be- 
tween the best-fit Ltr and that determined from 24 //m 
alone is similar to the 1 a uncertainties in Ltr. 

Using Spitze r and Herschel da ta for galaxies at red- 
shifts z < 2.0, lElbaz etlrtl (J2011D investigate the physi- 
cal reasons for the overestimation of Lir based on rest- 
frame Lg alone. These authors point out that the SED 
templates used to infer Lir are calibrated to match local 
ULIRGs, which are starbursts with compact, high sur- 
face density star formation. This contrasts with ULIRGs 
at high redshift (z > 2.0) which have more extended star 
formation occuring over longer timescales. These differ- 
ences in star formation surface densities and timescales 
lead to noticeable differences in the IR SEDs because 
of the variations in the spatial distribution of dust with 
respect to the massive stars that are heating this dust. 
In Section [5J we discuss the Lir/Ls ratios found here 
and place them in the context of the ratios found for 
other star-forming galaxies with similar luminosities but 
at lower redshifts. 

We conclude by noting that the iReddv et al.l (|2010f ) 
calibration, which is based on the correlation between 
L g and dust-corrected Ha luminosity at z ~ 2, pre- 
dicts Lir for LIRGs that are similar within the errors 
to those computed using the combined Spitzer, Herschel, 
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Fig. 5 . — C o mparison o f the ICharv fc Elbazl f2001). 

IDale fc Heioul J20021 ). IRieke et al] <200lD . and lElbaz et af] <20111) 
template fits to the observed 24, 100, and 160 fira and 1.4 GHz 
measurements, derived by scaling the template that best matches 
the observed fluxes, for L* galaxies at z ~ 2. The total infrared 
luminosities for the three templates are identical within the uncer- 
tainties (Table [4j. 

and VLA data. Finally, the infrared luminosities based 
on the radio data alone are in excellent agreement with 
the infrared luminosities inferred from fitting the dust 
templates to the observed 24, 100, and 160 /xm fluxes. 

4.4. Comparison of Dust SEDs 

Figure [5] compares the dust SEDs found for the 
luminosity-matched fitting to the 24, 100, and 160 /im, 
and 1.4 GHz data, weighted by their errors. The dust 
SEDs for the different models are broadly consistent 
with each other; one difference is in the stronger sil- 
icate absorption at 9/Ltm and additional aromatic fea- 
tures longward of 10 /im in the Rieke+09 model fit. The 
Rieke+09 templates are calibrated using Spitzer/IRS 
spectra, rather than photometry, between 5 and 36 /iin, 
and will understandably include mid-IR spectral features 
that are not present in the other models. 

The most notable difference between the templates can 
be seen in Figure [S] the Elbaz+11 and DH02 templates 
exhibit a broader range of dust temperatures with a 
colder dust component, relative to CEOl and Rieke+09. 
Resolving the full shape of the SED, and the average dust 
temperatures, will require larger stacked samples and 
deeper data in the submillimeter and millimeter regime. 
In any case, while some obvious differences remain be- 
tween the model templates given the data at our disposal, 
the integrals of these SEDs, namely the total Lir, are es- 
sentially identical within the uncertainties (Table 2]). For 
the subsequent analysis, we assume the Lir determined 
from the DH02 model that best fits the error-weighted 
fluxes at 1.4 GHz, 24, 100, and 160 /im. Assuming the 
Lir derived using any of the other models does not af- 
fect our conclusions. 

4.5. The Dust SED of Typical Star-Forming Galaxies at 

Z~2 

The average stellar population and dust SEDs for LJj V 
galaxies at z ~ 2 are shown in Figure [HI The optical pho- 
tometry indicates that an L^y galaxy at z ~ 2 has a UV 
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luminosity of Luv — 3.1 x 10 10 L© and the dust SED in- 
dicates a total infrared luminosity of Lm ~ 2.2 x 10 11 Lq. 
The Herschel and VLA data confirm the previous find- 
ing that UV-selected galaxies at z ~ 2 with 1Z < 25.5 are 
luminous in f rared galaxies (LIRGs; I Redd v et al.l 120051 
I2006U 120101 : lAdelberger fc Steid cl 2000); in particular, 
the medi an Ltr found here is virtually identical to that 
found bv lReddv et al.1 ()2010t ) based on an analysis of the 
24 /Am, Ha, and UV emission for UV-selected galaxies at 
z ~ 2. We also point out that the median value of Lir for 
LJjy galaxies is similar (to within a factor of rs 3) to the 
value of Lj R deduced f rom direct Spitzer de terminations 
of the IR LF at z ~ 2 (|Magnelli et al.ll2011[ ). 

The necessity of stacking the UV-selected galaxies, 
even for GOODS-dcpth Herschel PACS data, is illus- 
trated in Figure [7J Typical L* galaxies at z ~ 2 are eas- 
ily detected at wavelengths blueward of observed 24 /tin 
given the depths of the data considered here. Redward 
of this point, however, the dust emission is not sufficient 
for directly detecting these galaxies, thus stacking in the 
infrared and radio bands is required. Fortunately, the 
Herschel PACS sensivity and resolution are such that 
we can for the first time significantly detect the average 
thermal emission of non-lensed L* galaxies at z ~ 2. 

4.6. Variation of Lir with UV Slope and Bolometric 
Luminosity 

The stacks for the different subsamples indicate that 
galaxies with bluer UV spectral slopes are on average 
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cal (U n GTZ), ground-based near-IR (JK B ), Spitzer/IRAC (3.6 
8.0 /im), Spitzer/MlPS (24 /mi), Herschel/PACS (100, 160 /mi), 
and VLA 1.4 GHz data in the GOODS-North field, relative to the 
SED of an L* galaxy at z ~ 2. 

less infrared luminous than those with red UV spectral 
slopes (Table 0|. This systematic effect is approximately 
equal in magnitude to the uncertainties in the stacked 
flux measurements. Not surprisingly, those galaxies that 
are inferred to have iboi > 10 12 L based on the MIPS 
24 /im data (Sample E; Table 2]) have median stacked 
fluxes at 100 and 160 /im that are a factor of 3 — 4 times 
larger than the corresponding fluxes for the L* sample. 
The median infrared luminosity for Sample E is not as 
large as that computed from 24 /im alone and hence not 
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as large as the value of Lboi used to construct this sub- 
sample (e.g., Figure EJ. However, the best-fit template 
to the mid-IR, IR, and radio fluxes indicates that galax- 
ies in Sample E are still more infrared luminous than 
galaxies in other subsamples. These galaxies correspond 
to low-luminosity ULIRGs or higher luminosity LIRGs. 
In the next section, we compare the infrared and UV 
luminosities for each of the subsamples of UV-selected 
galaxies. 

5. DISCUSSION 

In the following, we discuss our results on L* galaxies 
at z ~ 2 in the context of the infrared properties of star- 
forming galaxies with similar luminosities at lower red- 
shifts. We then discuss the combined UV and IR lumi- 
nosity measurements for z ~ 2 galaxies and the implica- 
tion for their average dust attenuation. We also compare 
the measured dust attenuation with that inferred from 
the local correlation between dustiness and UV slope. 
Finally, we use the Herschel data to determine the me- 
dian bolometric luminosities and star formation rates of 
galaxies in our sample, the correlation between bolomet- 
ric luminosity and dust attenuation, and the variation in 
dust attenuation with UV luminosity. 

5.1. Ratio of Infrared to Mid-Infrared Luminosity 

For a consistent comparison with lElbaz et al.l (|2011h . 
we have recomputed Lg using the mid-IR SED of M82 to 
fc-correct the 24 /an flux. The implied fc-corrections are 
not substantially different than those obtained from the 
average mid-IR SED of the 1 2 local star-forming galaxies 
listed in lReddv etall (|2006bf ). Adopting the CE01 lumi- 
nosity matched and weighted value of Lir, we compute 
IR8= 7.7±1.6, 8.9±1.3, 8.3±2.4, 8.9±1.4, and 8.4±0.5, 
for Samples A, B, C, D, and E, respectively. Assuming 
the Lir computed from the 100 and 160 /im data only re- 
sults in ratios that are not significantly different than the 
ones quoted above. These values imply that the galaxies 
in our sample predominantly lie on the boundary be- 
tween the ratios found for infrared main sequence galax- 
ies (IR8= 4.9^22)) an d those found for starburst galaxies 
which have an exponential tail distribution extending to 
IR8= 15 — 20. The relatively high IR8 ratios also may 
explain why the typically-adopted templates (e.g., CE01) 
tend to not fail as badly for galaxies in our sample rela- 
tive to more luminous ULIRGs at z y 2 when inf erring 
L IR from L 8 alone (e.g- lReddv et al.ll2006bll2010D . 

A detailed comparison between the morphologies of the 
z ~ 2 galaxies studied here and local star-forming galax- 
ies is beyond the scope of this paper. However, we can 
still make some general inferences regarding the degree of 
"compactness" of the IR emission in these galaxies and 
how it may affect their IR8 ratios. We compute the I R 
luminosity surface density following lElbaz et al.l (|2011l ): 



Sir = 



Ljr/2 



(1) 



'IR 



The typical UV half-light radius of galaxies in our sam- 
ple is njv — 2kpc. If we assume that the IR half- 
light radius is rou ghly half this value, as suggested by 
lElbaz et"aL1 ([201 ID . then r m « lkpc, implying Sir w 
3 x lO lo L0kpc~ 2 . On the other hand, if hr w njv 
as might be expected for galaxies at high redshift where 



the UV emission is dominated by OB stars, then Sir « 
8 x 10 9 L kpc -2 . Even this lower value is a factor of 4 
larger than the Sir < 2x lO 9 L kpc~ 2 typical of the 
infrared main sequence ga l axies with extended star for- 
mation in the lElbaz et al.1 (|2011l ) sample. The relevance 
for the present discussion is that while L* galaxies at 
z ~ 2.3 appear to have IR8 ratios that are a factor of sa 2 
larger than for main sequence galaxies at lower redshifts, 
they also exhibit Sir that are larger than those found 
for main sequence galaxies, and hence are more compact 
for their IR luminosity. Formally, galaxies in our sample 
would lie on the "boundary" between IR main sequence 
and IR starburst galaxies. These larger IR8 ratios result 
in the over-prediction of Lir based on Lg alone when 
using the standard templates, as noted in Section [4] 

One possibility is that the IR8 ratio for typical star- 
forming (L*) galaxies may evolve with redshift, transi- 
tioning from IR8 ~ 4.9 at z < 2.0 to IR8 ~ 8 at z ~ 2.3. 
This effect is likely due to the smaller sizes of high red- 
shift galaxies for their IR luminosities relative to lower 
redshift galaxies. Finally, it is worth noting that the 
star formation rate surface density appears to be roughly 
constant for L* galaxies at high er redshifts (z ~ 4 — 7; 
S SFR ~ l^Moyr^kpc- 2 ; e.g.. lOesch et alj|2010t ). sug- 
gesting that there may be a similar commonality in the 
IR8 ratios (if one could measure them) for typical star- 
forming galaxies at very high redshift, just as is observe d 
for main sequence galaxies at z < 2.0 (El baz et al.ll201ll) . 

5.2. Dust Obscuration of Typical Star-Forming Galaxies 

at z ~ 2 

5.2.1. Definitions Relevant to Dust Obscuration 

Before proceeding, it is useful to define several terms 
that have been typically used interchangeably in the 
literature. First, we define "dust obscuration", or at- 
tenuation, as Lir/Luv- Note that this ratio is not 
equivalent to the ratio of obscured to unobscured star 
formation rate, SFRir/SFRuv, given the difference in 
scaling required to convert the UV and IR luminosi- 
ties to star formation rates. We also define the "dust 
correction factor" needed to recover the total star for- 
mation rate from that computed based on the unob- 
scured UV luminosity as (SFRir+SFR uv )/SFR uv = 1+ 
SFRir/SFRuv- For the L* sample, the median dust ob- 
scuration is Lir/Luv — 7. 1 ± 1 . 1 , the ratio of obscured to 
unobscured SFR is SFRir/SFR uv = 4.2 ± 0.6, and the 
dust correction factor is 5.2 ± 0.6. The L\jy, Lir, ratio 
of obscured to unobscured SFR, and total SFR for each 
subsample are list ed in Tab l e \5\ For the conversion to 
SFR, we assume a lSalpete"n (|1955l ) ini tial mass function 
with limits from 0.1 to 100 M and the lKennicuttl (|1998l ) 
conversions betweeen UV/IR luminosity and SFR. The 
results indicate that roughly 80% of the star formation 
is obscured for L* galaxies at z ~ 2. 

The dust obscuration varies between Lir/Luv — 4.8 ± 
1.0 for the subsample with blue UV slopes (Sample C) 
and Lir/Luv = 10.7±1.9 for the subsample with red UV 
slopes (Sample D), and is as high as Lir/Luv ~ 37.7±7.9 
for the most bolometrically-luminous galaxies (Sample 
E). These observations imply a trend in dustiness with 
both UV slope and bolometric luminosity. 

5.2.2. Validity of the UV Attenuation Curve for L* Galaxies 

at z ~ 2 
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TABLE 5 
Properties of the Stacks IV: Derived Quantities 





iuv a 


L m b 






SFR(UV)±SFR(IR) d 


Sample 


(10 n L Q ) 


(lO 11 ^) 


1±SFR(IR)/SFR(UV) 


cd 


(Meyr 1 ) 


A. 


0.32 ±0.02 


2.3 ±0.3 




5.3 ±0.6 




49 ±6 


B. 


0.31 ±0.02 


2.2 ±0.3 




5.2 ±0.6 




47 ±6 


C. 


0.33 ±0.03 


1.6 ±0.3 




3.9 ±0.6 




37 ±6 


D. 


0.28 ± 0.02 


3.0 ±0.5 




7.2 ±1.1 




60 ±9 


E. 


0.35 ±0.07 


13.2 ±0.8 




20.9 ±3.2 




238 ± 15.5 


F. 


0.38 ±0.05 


3cr: < 2.0 




3cr: < 2.4 




3cr: < 58 


a Mean and error in mean 


of UV luminosity 


in units of 10 Lq. 






Infrared 


luminosity, in 


units of 10 L0, derived from color-matching 
els to the observed fluxes. For Sample F, we as 


and normalizing the 


IDale & Heloul (120021) mod 


iume the upper limit 



in Lift implied by the observed fluxes at 24 and 160 //m and the upper limit at 100 /im and 

1.4GHz. 

c Dust correction factor required to recover the total SFR from the UV-dctcrmincd SFR. 

Median star formation rates in Mg yr _ assuming a Salpctcr (1955) IMF from 0.1 to 100 Mq 
and the Kcnnicutt (1998) relations between UV/IR luminosity and star formation rate. For the 
"young" subsamplc (Sample F), we multiply the UV SFR determined from the Kcnnicutt ( 1998) 
relation by a factor of 2. This is done to account for the fact that the mix of O and B stars 
contributing to the UV continuum emission has not equilibrated for ages < 100 Myr (assuming 
a constant star formation); thus, the Kcnnicutt (1998) conversion between UV luminosity and 
SFR will undcrprcdict the total SFR for such "young" galaxies. 



In particular, we show the dust obscuration derived 
for these samples as a function of UV slope, j3, in Fig- 
ure [3 Up to luminosities of Lm w 1O 12 L , we find 
that galaxies with redder (3 are on average dustier. Fur- 
thermore, the correlation between dustiness and UV 
slope is essenti ally identical to that found for local star - 
burst galaxies (jMeurer et al.lll999t ICalzetti et al.ll2000D . 
This result has been found by several other inves- 
tigations targeting moderately luminous galaxies and 
using a variety of star formation tracers at z ~ 2 
(e.g., IReddy fe Steidel 12 004: R eddvet all l2006bL l20l(it 
Daddi et al. 2007: Pannclla et al. 20 091) a n d z ~ 3 (e.g. 
Seibert et all 120021 : iNandra et al.l 120021 : iMagdis et ail 
2010aHbf ). 



Infrared selection, e.g., such as 24 ^m selection, gener- 
ally results in samples where the b ulk of galaxies do not 
abide by the lMeurer et al.l(ll999D orlCalzetti et al.l (|2000l ) 
attenuati on curves (e. g ., IMurphv et al.l [20111). As dis - 
cussed in lReddvetail (|2006bl ) and IReddy et al.l (|2010D . 
the correlation between UV slope and dust attenua- 
tion breaks down for more infrared luminous galaxies at 
z ~ 2, as well as younger galaxies with ages < 100 Myr 
at the same redshifts, where the latter tend to follow a 
steeper attenuation curve (i.e., they ar e less reddened at a 
given UV slope than predicted by the lMeurer et al.lll999l 
relation). Systematic deviations from the local star burst 
attenuation curve have also been observed at lower red- 
shift (z < 2) based on Herschel/PACS and S PIRE data 
(e.g.. lBuat et al.l|2010t iBurgarella et al.ll2011l) and Akari 
data (|Buat et al.ll2011l) . 

Much of the aforementioned deviation from the local 
starburst attenuation curve can be understood in the 
context of the range of bolometric luminosity probed by 
the different UV and IR selections. UV color selection is 
sensitive to galaxies with moderate [L* ) luminosities and 
lower dust extinction than those selected in the infrared. 
Because dust attenuation is a strong fun ction of bolo- 
metric luminosity and the validity of the [ Meurer et al.| 
(Il999f) relation is luminosit y dependent (Meurer et al. 
fl999l iGoTdader et all l2003 iReddv et~aTl l2006bl . I2010Q . 
it is natural to expect departures from this relation for 
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Fig. 8. — Mean dust attenuation (Lir/Luv) versus UV slope 
(/3) for different subsamples of z ~ 2 galaxies. Also shown are 
attenuation curves f or the SMC and for local UV starbursts from 
[Meurer ct al. (1993), an d t ne 3<r upper limit and stacked 24 (im 
implied value {cyan point) of the dust attenuation for the youngest 
galaxies in our sample. 

galaxies that may be selected via their infrared emission. 
Deviations may also be observed in infrared luminous 
galaxies that also have large stellar masses; in this case, 
the UV continuum associated with the massive OB stars 
may be extinguished relative the UV emission from less 
massive stars, resulting in a redder UV continuum slope 
for a given dust obscuration (e.g.. IMurphv et all l2011t 
iBuat et aLll2010D . These results stress that one must 
take care in applying any starburst attenuation curve 
to high redshift galaxies, with the realization that such 
relations may apply in one regime but fail in another de- 
pending on the properties of the galaxies in one's sample. 
Here, we have shown that galaxies that lie around L* of 
the UV luminosity function have dust obscuration - as 
measured from Spitzer, Herschel, and VLA data - that 
correlate with their UV slopes, and that this correlation 
is similar to that observed for local starburst galaxies. 



14 



The stacked Herschel data do not directly indicate the 
intrinsic dispersion in the relation between UV slope and 
dustiness. However, an indirect estimate of this scatter 
comes from an analysis of the 24 /im data. Specifically, 
109 of 311 galaxies i n the larger (and multiple field) sam- 
ple of iReddv et al.l (J2010I ) are detected at 24 /zm. Using 
a survival analysis t o take into a c count both detections 
and non-detections, IReddv et al.l ()2010D found a disper- 
sion of w 0.40 dex between dust attenuation, L m /L V y, 
and rest-UV slope, j3. If the scatter in th e Ls-to -LrR ratio 
is similar to that found by lElbaz et al.l (|2011h for lower 
rcdshift (z < 2.0) galaxies with Lir < 10 12 L (having a 
la dispersion of « 0.1 dex), then the implied total dis- 
persion in the relation between dustiness and UV slope is 
f* 0.45 dex. The correspondence between Lir/Luv an d j3 
to within a factor of » 3 thu s verifies the applicab ility of 
the lMeurer et al.l (|1999h and lCalzetti et al.l <|2000ft atten- 
uation curves for galaxies with Lir < 10 12 L at z ~ 2. 
The UV attenuation curve is sensitive primarily to the 
geometry of dust and stars within galaxies, and/or vari- 
ations in dust composition. Hence, the correspondence 
of the UV attenuation curves between local starbursts 
and z~2L* galaxies implies a remarkable similarily in 
the processes that give rise to the spatial distribution of 
dust and stars and the dust composition in galaxies over 
w 10 billion years of cosmic history. 

5.3. Comparison with Recent Studies 

Using the Herschel data to directly probe the thermal 
dust emission, we have shown that the local correlation 
between UV slope and dust obscuration remains valid for 
typical star-forming galaxies at z ~ 2. In the following, 
we compare our results with several recent studies of the 
dust attenuation of high redshift galaxies. 

5.3.1. Radio Emission from UV-selected Galaxies 

ICarilli et al.l (|2008f) stack the radio emission of z ~ 
3 LBGs in the COSMOS field and deduce that 
SFR radio /SFR uv = 1.8 ± 0.4, wh ere SFR radio is derived 
using the calibration of lYun et~a l. (2001). This calibra- 
tion is based on equating the local star formation rate 
density with the integral of the radio luminosity function, 
so in this case the radio SFR should represent a total SFR 
(but see below) , including obscured and unobscurcd com- 
ponents. The factor of 1.8 from the ICarilli et al.l (2008) 
study is significantly smaller than the factor of ~ 5 com- 
p uted for our UV-sel ected sample at z ~ 2. 

ICarilli et al.l (|2008l ) discuss several possibilities for the 
suppression of radio flux with respect to SFR at z ~ 3. 
For a consistent comparison, we assess their results us- 
ing the same calibratio n used in t his analysis. In par- 
ticular, we employ the iBelll (|2003f) correlation between 
total infrared luminosity (Lir) and specific luminosity 
at 1.4 GHz. Doing so, the stacked median radio lumi- 
nosity of LBGs in the COSMOS field, L 1A = 5.1 x 
10 29 ergs -1 Hz -1 , corresponds to Lir rs 2.2 x 10 11 L Q . 
These values are essentially identical to those determined 
for our L* sample at z ~ 2: L1.4 = (5.2 ± 1.0) x 
10 29 ergs" 1 Hz" 1 and L m = (2.2 ± 0.3) x 10 11 L Q (Ta- 
ble HI). The obscured SFR correspo nding to the Lir 
for the COSMOS LBGs, assuming the IKennicuttl (|l998l ) 
relation, is w 38M yr _1 . Hence, the factor that we 
compute to recover the total star formation rate from 



the UV star formation for the ICarilli et al.l (|2008f ) sam- 
ple, assuming their value of the unobscured SFR of 
lTMoyr" 1 , is 1 + 38/17 « 3.2. This is clos e to a factor 
of two larger than the value of 1.8 given in ICarilli et al.l 
( 2008 ) . This d i screp ancy results from the fact that while 
the lYun et all (|2001l ) correlation between radio luminos- 
ity and star formation rate gives a star formation rate 
that is in good agreement with that inferred from Lir 
(SFR~ SSMgyr" 1 ), it underestimates the total SFR 
(SFRir+SFRuv) of ~ 55 M Q yr" 1 . 

To illustrate this point, we plot in Figure [9] the rela- 
tionship between total st ar formatio n rate and infrared 
luminosity, adopting the IKennicuttl (|1998| ) conversions 
between UV/IR luminosity and sta r formation rate, fo r 
the sample of 392 z ~ 2 ga laxies of | Reddv et a l. (2010), 
and for the local samples of Bell ( 2003) and Huang et al. 
(2009). The top axis in Figure [3] indicates the radio lumi- 
nosity that corresponds to Lir assuming the lYun et al.l 
(2001) calibration, and the solid line shows the relation- 
ship between radio lumino sity and star for mation rate 
derived in that study. The lYun et al.l (|2001l ) calibration 
is valid if most of the star formation is obscured, as is 
the case for ULIRGs at both z ~ 2 and z ~ 0. However, 
a substantial fraction of the LIRGs in the z ~ 2 sam- 
ple have a significant contribution from unobscured star 
formation, where the unobscured component is at least 
50% of the obscured star formatio n. The same is tru e 
for the COSMOS LBG sample of ICarilli et all (120081) . 
The ratio of obscured to unobscured star formation rate 
is a strong function of total star formation r ate or bolo- 
metric luminosity (e.g., Figure^ Rc ddv et al.l l2010). and 
the UV component obviously cannot be neglected for ob- 
jects that have significant UV emission. Figure [5] shows 
that one would significantly underestimate the total star 
formation rate of LIRGs at z ~ 2 based on their IR emis- 
sion alone (the lKenn icutt 1998 relation between SFR and 
Lir is only valid in the opti cally thick limit), or based 
on using the lYun et al.l (|2001l) relationship between radio 
luminosity and star formation rate. 

Note also that Carilli et al. (2008) combine the median 
radio SFR with the mean unobscured UV SFR to deter- 
mine the effect of dust. In general, mean luminosities will 
b e larger th a n med ian ones for a population drawn from 
a iSchechteH (|1976l) luminosity function, modulo sample 
incompleteness. For our L* sample, the mean UV lu- 
minosity is about 15% larger than the median. More 
importantly, the m ean unobscured UV luminosity of the 
ICarilli et all (|2008f ) sample of 17 M Q yr _1 is about a fac- 
tor of two larger than the median UV SFR for our L* 
sample of rs 8.5M yr _1 . I t is possible that t he rela- 
tively UV-bright LBGs of the lCarimlTaTI (|200l sample 
are somewhat less attenuated than more typical (and less 
UV-luminous) LBGs at z ~ 2 — 3, though this is contrary 
to what h as been found for U V-selected galaxies at these 
redshifts (Red dv et al.ll20l6T ). Without further analysis 
of their candidates, we c onclude that the low er UV obscu- 
ration factor deduced bv lCarilli et al.l (2008) is likely due 
to the brighter unobscured UV luminosity of their can- 
didates, relative to the obscured star formation. What is 
clear from our spectroscopic sample is that the stacked 
radio flux for UV-selected galaxies at z ~ 2 predicts an 
Lir that is identical to that obtained using direct mea- 
surements of thermal dust emission from the Herschel 
data, and this value of Lir implies a dust correction of a 
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Fig. 9. — Total star formation rate versus infrared luminos- 
ity (Lm) f° r the sample of 392 UV-selected z ~ 2 galaxies of 
I Redd v ct al. (20 13) ( r ed cir cles a nd orange up p er lim its), and the 
local samples of Bell (2003) and Huang ct al. (2009) (dark green 
circles). The top axis shows the radio lumi nosity that corres ponds 
to Ljr assuming the radio-IR correlation of Yuri ct al. (2001), and 
the solid line denotes the relationship between radio luminosity and 
star formation rate derived in that study. The shaded region covers 
the area where the unobscured star formation is at least 50% of the 
obscured star formation. Th e large open s tar de notes the position 
of the COSMOS LBGs from lCarilli et al.l 12003 1. 

factor of « 5. 

5.3.2. Investigations of the Extragalactic Background Light 

(EBL) 

A second and more recent study that has suggested 
obscuration fa ctors that are differ ent from those pre- 
dicted from the iMeurer et al.l (|1999f ) relation comes from 
an analysis of the infrared backg round based on s tacke d 
Spitzer and BLAST data by ICharv fc Popl (|2010f ). 
These authors estimate the extragalactic background 
light (EBL) contributed from galaxies at z > 1 and 
conclude that the UV-based dust corrections for typical 
star- forming galaxies (e.g., those selected by th e ir UV 
emission) must be lower than the IMeurer et al.l (|1999f ) 
prediction in order that their integrated emission does 
not violate the background light constraints. However, 
directly probing the thermal dust emission of L* galaxies 
at z ~ 2, as we have done here, shows conclusively that 
these galaxies have dust attenuations tha t are similar, 
on ave rage , to those computed b ased on the lMeurer et al.l 
(1999) and lCalzetti et all (|2000D attenuation curves (Fig- 
ure |5J). How can we reconcile these two results? 

The EBL is sensitive to the average dust attenuation 
of all g alaxies, not just the UV- brig ht ones studied here . 
Indeed. Ileddv fc Steidell rf2009h and lReddv et~all (|2"o"l0h 
use physical arguments and stacked Spitzer data to show 
that the dust obscuration of UV-faint galaxies is lower 
than in UV-bright ones (see also next section), and that 
this luminosity dependence implies that the average dust 
obscuration, integrated over the entire luminosity func- 
tion, can be close to a factor of two lower than the mean 
dust obscuratio n found for just the UV bright galaxies 
(see Table 5 of iReddv fc Steidell 120091) . This does not 
necessarily imply a failure of the IMeurer et al.l (|1999| ) re- 



lation for typical star-forming galaxies at high redshift; it 
simply means that the rest-frame UV slope also becomes 
bluer, on average, for UV-faint galaxies, as the empiri- 
cal ev idence seems to indicate at z ~ 2.5 (|Bouwens et al.1 
2009). What is clear from the present analysis is that 
the direct measurements of the thermal dust emission of 
UV-selected galaxies at z ~ 2 imply that the local UV 
attenuation curve remains valid for these galaxies. 

5.3.3. Stacked ¥Lerschel/SPIRE Measurements for z ~ 2 
UV-selected Galaxies 

Strong source confusion in the Herschel/ 'SPIRE 250, 
350 and 500 /im data make it very difficult to carry out 
stacking experiments capable of reaching flux limits as 
faint as those we expect for L* UV-selected galaxies at 
z ~ 2 (approximately 1 mJy - see Figure |5]). After some 
experimentation, we have chosen not to include those 
data in our analysis. Nevertheless, iRigopoulou et al.l 
(2010) stack SPIRE 250 /im data for a smaller sam- 
ple of brighter, 24 /xm-detected, UV-selected galaxies at 
z ~ 2. These were taken from the sample of lReddv et al.l 
(|2006al) , also used here, but were limited to a 69 ob- 
jects with individual 24 /mi detections. That subsamplc 
is therefore likely to be more IR-luminous, on average, 
than the larger sample considered he re, where individ- 
ual de tection at 24 /jm is not required. IRigopoulou et ahl 
(2010) measure a stacked 250 /im flux of /250 = 2.7 ± 
0.8 mJy, corresponding to a total infrared luminosity 
L m « 4.2 x 10 11 L at (z) « 2. 21 This value is sim- 
ilar to what we predict at 250 /im when we stack the 
100 and 160 /im data for the same set of galaxies (i.e., 
the brighter objects with individual 24 /im detections). 
Given that we have (1) controlled for issues of clustering 
and confusion at 100 /im, and cross-checked the 160 /im 
fluxes with those obtained at 100 /im, and (2) performed 
detailed simulations to demonstrate the robustness of 
the stacked fluxes (Section [3]), we believe the PACS con- 
straints on the median IR SED to be robust. With the 
PACS data we are able to constrain the median IR SED 
for L* galaxies with lower average infrared luminosities, 
irrespective of whether they are individually detected at 
24 /im. Future observations with the Atacama Large Mil- 
limeter Array (ALMA) should provide higher resolution 
and more sensitive observations of the Rayleigh- Jeans 
emission from typical star-forming galaxies at high red- 
shift. 

5.4. Dependence of Dust Attenuation on Stellar 
Population Age and Luminosity 

5.4.1. Galaxies with "Young" (< 100 Myr) Stellar 
Population Ages 

Roughly 14% of the galaxies in our entire sample are 
identified as having a "young" stellar population with in- 
ferred ages of < 100 Myr. Stacking these galaxies results 
in a formal nondetection (S/N < 3) in the Herschel and 
VLA stacks (Table[2]). The implied 3 a upper limit to the 
dust obscuration suggests that these galaxies follow an 

21 We note that the stacked 250 (im flux shown in Figure 2 of 
IRigopou lou ct al. (2010) appears to be nearly lOx fainter than the 
value cited in the text. The authors state that they derive a total 
infrared luminosity Lib = (1.5±0.5) X 10 Lq using the CE01 dust 
templates, but we are unable to reproduce this, finding instead a 
luminosity that is 2.8 X larger, based on the value /250 = 2.7 mJy 
cited in the text. 
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Fi g. 10. — Left: Bolo metric luminosity, Lbol = LlR + Ltjv, as a function of Ltr/-Ltjv for a sample of local galaxies from Bell (2003) 
and Huang ct al. (2009), shown by the open diamonds. The sample of 392 UV-selected galaxies at 1.5 < z < 2.6 analyzed in Re ddv et al.l 
(120101 ) are shown by the open circles and upper limits. The relationship between Lbol and Lir/Ltjv for Samples A, C, D, and E, based 
on the stacked Herschel, Spitzer, and VLA data, are denoted by the large filled squares (error bars are not shown for clarity). The same 
quantity for the L* sample at z ~ 2 (Sample B) is show n by the large open square. Right: Correlation between Lbol an d Ltjv, where the 
cyan points denote the results using 24 /im data and the[Rcddy et al. (2010) calibration between Lg and Ljj^. The purple points indicate 
results from incorporating the Spitzer/MIPS 24 fim, Herschel/PACS 100 and 160 /im, and VLA 1.4GHz data to determine the median IR 
SED as a function of UV luminosity. 

attenu ation curve that is " steep er" than the typically as- 
sumed [MeureLeES] (| 19991) and lCalzetti et all (|2000h at- 
tenuation curves (Figure[5]), in the sense that they exhibit 
a redder UV slope for a given obscuration. More strin- 
gent constraints on the dust attenuation of these galaxies 
conies from their stacked 24 /jm emission. Assuming the 
DH02 template that best matches the median 24 /im flux 
of the "young" galaxies implies Ltr/£uv ~ 1-55 ± 0.35. 
The median UV slope of these galaxies of j3 = —1.30 
then implies that their dust obscurations are consistent 
with the SMC attenuation curve. The deviation of such 
young, high-redshift galaxies from the standard a ttenua - 
tion curve is a result first noted in Rcddy ct al. (2006b) 
and further investigated in iReddv et al. (2010). IRS 
spectrocopy of a couple of young lensed LBGs provide ad- 
ditional evidence tha t such young galaxies ar e less dusty 
at a given UV slope (jSiana et al.ll2008l |2009|) than their 
older and m ore typical counterparts at the same redshifts 
(jReddvet alj[20ia) . 

The small carbonaceous dust grains giving rise to the 
mid-IR emission are believed to be produced in AGB 
stars over a longer timescale than the Type II SNe re- 
sponsible for the large grains that emit thermally in the 
infrared (e.g., Gall iano et al.l [20081) . It is therefore rele- 
vant to ask whether evolution in the small-to-large dust 
grain ratio may be responsible for the observed differ- 
ences in attenuation curve with stellar population age. 
Based on an X-ray stacking analysis and an examina- 
tion of_the_theOTeticaU PAH metallicity from dust mod- 
els, IReddv et al.l (|2010l ) argued that the small-to-large 
dust grain ratio is unlikely to change significantly over 
the relatively narrow range in metallicity probed by the 
UV-selected sample. IRS spectroscopy of at least a cou- 
ple of young lensed LBGs at high-redshift also indicates 
that their ratios of mid-to-far IR emission are similar 
to those found for local starburst galaxies (jSiana et al.1 



120081 120091) and similar to the rat ios observed for z ~ 2 
galaxies (|Reddv et al.ll2"006bLl2010t ). The Herschel obser- 
vations, which probe the thermal dust emission at z ~ 2, 
are consistent with these findings based on 24 /im emis- 
sion alone. 

If the difference in attenuation curve is due purely to 
geometrical effects, then it suggests that the dust cover- 
ing fraction of high-redshift galaxies evolves significantly 
with time. A possible simple scenario is one in which the 
first generation of stars quickly pollutes the ISM with 
dust and metals over a dynamical timescale of a few tens 
of millions of years, at which point much of the dust is 
foreground to the stars and the dust covering fraction 
is high. As star formation proceeds and either becomes 
spatially extended or is able to drive outflows sufficient 
in momentum to perturb the dust/gas distribution of the 
ISM, the resulting attenuation becomes patchier. So, we 
might expect in this situation that galaxies would grad- 
ually transition from a steep attenuation curve like that 
of the SMC to a grayer starburst attenu ation curve (e.g. 



iMeurer et aT]|199£ 



Calzetti et al.ll2000l see also discus- 
sion in iBuat et al.l 1201 ll ) . Detailed studies of the vari- 
ation of the interstellar absorption lines (as a proxy for 
the dust covering fraction) with dust attenuation will be 
needed to test this scenario. 

5.4.2. Variation in Dust Attenuation with Bolometric and 
UV Luminosity 

MIPS 24 /an studies have shown that z ~ 2 galax- 
ies follow a tight trend be tween bolometric luminosity 
and dust attenuation (e.g., IReddv et al.ll201CJI and refer- 
ences th erein), similar to that observed locally a nd at 
z ~ 1 (iBurgarella et all l2009t IBuat et al.l [20071 120091: 
IWang fc Heckmanl Il996f ). However, the normalization 
of this trend depends on redshift, such that at a given 
bolometric luminosity, galaxies at high redshift are on 
average less dusty than local ones (|Reddv et al.l l2006bl 
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l2008ll2010D . This result is not unique to our UV-selected 
sample. Rest-frame optical, near-IR, and submillimeter 
selected galaxies at z ~ 2 also show a similar offset in 
obscuration per unit star format io n compared with lo - 
cal galaxies (| Redd vet all l2006bl ) . IReddv et all (|2010l ) 
demonstrate that this relationship is likely driven by 
metallicity evolution. The stacked Herschel and radio 
data confirm these previous results (Figure ITUl) . This is 
not surprising given the similarity in Lir and dust at- 
tenuation inferred in this study with those inferred for 
L* galaxies based on 24 /jm data alone. The combined 
Herschel, Spitzer, and VLA data confirm that £{j V S a l ax ~ 
ies at z ~ 2 are a factor of » 10 less dusty than galaxies 
with similar bolometric luminosities in the local universe. 
Similarly, L{j V galaxies at z ~ 2 are a factor of « 20 — 30 
times more bolometrically-luminous than local galaxies 
with a similar dust obscuration (Figure [TU]) . These evo- 
lutionary effects can also be seen in Figure [9] In the 
local universe, it is only for galaxies with Lir < 10 11 L Q 
(i.e., galaxies fainter than LIRGs) where the unobscured 
star formation begins to contribute significantly to the 
total star formation rate. In contrast, most lower lu- 
minosity LIRGs at z ~ 2 have a significant fraction of 
unobscured star formation, as indicated by the plume of 
z ~ 2 galaxies extending away from the SFRboi =SFRir 
line in Figure HO indicating that LIRGs at z ~ 2 are more 
UV-transparent relativ e to LIRGs in the lo cal Universe. 
Finally, we note that I Redd v et al.l (|2010() use a 24 /im 
analysis of a larger sample of UV-selected galaxies at z ~ 
2 to show that the fraction of 24 /im detections of these 
galaxies decreases by a factor of two proceeding from UV- 
bright (L uv « 1O U L ) to UV-faint (Luv » 1O 1O L ) 
galaxies (see Figure 14 in IReddv et~aLll2010D . This UV 
luminosity trend in 24 /im detection fraction implies that 
UV-faint galaxies are on average less infrared-luminous 
than their UV-bright counterparts, a result that is fur- 
ther confirmed by stacking the 24 /im emission of galax- 
ies in bins of UV luminosity. A stacking of the combined 
Spitzer, Herschel and VLA data as a function of UV lu- 
minosity shows a similar trend, indicating that UV-faint 
galaxies are also less bolometrically luminous than UV- 
bright galaxies ( Figure [TO]). The actual dust attenuation, 
Lir/Luv, is roughly constant (within the uncertainties) 
with decreasing UV luminosity, as both UV and IR lu- 
minosities decrease in tandem. In any case, it is clear 
from these observations that care must be taken when in- 
ferring the average dust obscuration and Lir, and their 
effect on global quantities like the star formation rate 
density, gi ven the UV luminosity dependence of th ese 
quantities (|Reddv et alJl2008t IReddv fc Steide]|l2009t ). 

6. CONCLUSIONS 

We have used the deep 100 and 160 /jm data of the 
GOODS- Herschel Open Time Key Program, supple- 
mented with deep Spitzer/MIPS 24 /im and VLA 1.4 GHz 
radio imaging, to investigate the infrared luminosities 
and dust obscuration of typical star-forming galaxies at 
high redshift. We focus on the median stacked mid- 
infrared, far-infrared, and radio fluxes of a sample of 
146 UV-selected galaxies with spectroscopic redshifts 
1.5 < z spcc < 2.6 in the GOODS-North field. These 
galaxies have luminosities around L* of the UV luminos- 
ity function at these redshifts. 

Because these galaxies are individually undetected at 



100, 160 /im, and 1.4 GHz, we perform median stacking 
analyses to measure their average fluxes. Three tests are 
performed to verify the robustness of our stacked results. 
The first test measured the effects of clustering. The sec- 
ond test determined the chance probability of measuring 
a stacked flux as high as the one obtained for the target 
galaxies. The third test allows us to measure the bi- 
ases and uncertainties in stacked flux measurements by 
stacking on the positions of artificial sources added to 
the images. 

To interprete these fluxes, we consi der a vari- 
ety of dust templates , including those of Elbaz et al.1 
d20Tl.lCharv fc ElbadlpOOlh . lDale fc Heloul l|2002) ■ and 
iRieke et al.l (|2009l ) , as w ell as the 24 /im-Lm and r adio- 
infrare d correlations of IReddv et al.l ((2010) and iBelll 
(2003), respectively. Fitting these templates to the bias- 
corrected fluxes and infrared colors reveals that LJj V 
galaxies at z ~ 2 with UV luminosities Luv ^ 10 10 Lq 
have a median infrared luminosity of Lir = (2.2 ± 0.3) x 
10 n L Q . Galaxies in our sample exhibit Lir/L§ (IR8) 
ratios that are a factor of rs 2 larger than those found for 
most star-forming galaxies at lower redshifts, likely due 
to the fact these UV-selected galaxies are relatively com- 
pact for their infrared luminosity. One possibility is that 
the roughly constant IR8 ratio observed for most (main 
sequence) galaxies with redshifts z < 2.0 likely shifts to- 
wards larger values at z > 2.0, due to the fact that the 
galaxies are on average smaller for their IR luminosity 
relative to lower redshift galaxies. 

Stacking the Spitzer, Herschel, and VLA data as a 
function of UV spectral slope, /3, and bolometric lumi- 
nosity, Lbolj indicates that galaxies with redder /3 and 
higher Lboi have on average larger infrared luminosities, 
in accord with expectations based on trends found for 
local star-forming galaxies. Based on the Lir, we pro- 
ceed to examine the dust attenuation, Lir/Luv > and 
the dust correction factor (i.e., the factor required to 
recover the bolometric SFR from the unobscured UV 
SFR), 1+ SFRir/SFRuv, for galaxies in our sample. 
For L* galaxies at z ~ 2, we find a dust attenuation 
of Lir/Luv = 7.1 ± 1.1, which corresponds to a dust 
correction factor of 5.2 ± 0.6, implying that ~ 80% of 
the star formation is obscured. This result is consis- 
tent with those found in previous studies of the dust- 
corrected UV and Ha, 24 /jm, and radio and X-ray stacks 
of the same galaxies examined here (jReddv fc Steidell 
[200llReddv et al.ll2006bt l2010h . 

We have examined the relationship between the UV 
spectral slope, /3, and dustiness for z ~ 2 galaxies. A 
small fraction (w 14%) of galaxies are identified as hav- 
ing stellar population ages < 100 Myr from fitting stellar 
population SEDs to the observed optical through near- 
IR photometry. The upper limit in dust attenuation for 
these "young" galaxies suggests that they may follow a 
"steeper" attenuation curve than the one obse rved for lo- 
cal starburst galaxies (jMeurer et al. 1999; Calzctt i et al.l 
2000), in the sense that they are less dusty (i.e., have 
lower Lir/Luv ratios) for a given /3 than older and more 
typical galaxies at the same reds hift. This result was 
first noted in the 24 /im an alyses of Rcddy ct al. (2006b) 
and further investigated in IReddv et al.l (|2010D . and may 
be attributable to the higher dust covering fractions in 
young galaxies. When considering more typical galaxies 
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with ages > 100 Myr, we find that their median -Lir/I/tjv 
increases as /3 becomes redder, and that this correlation 
is essentially identical to that found for local starbursts 
(|Meurer et 'aIJI1999tlCaIzetti et al.ll2000j ). Comparison of 
the Herschel results at z ~ 2 with measurements of local 
galaxies confirms the previously found trends that im- 
ply that LIRGs at z ~ 2 are more UV transparent (i.e., 
less dusty) than LIRGs in the local universe, and that 
UV-faint galaxies at z ~ 2 have lower Lm and hence 
fainter bolometric lumin osities than UV-br ight galaxies 
at z ~ 2 (e.g., see also iReddv et al.l 120101 ) . Based on 
the direct Herschel measurements of the rest-frame ~ 30 
and 50 /im thermal dust emission, we find that the local 
UV attenuation curve holds at z ~ 2 for galaxies with 
Lboi ^ 10 12 L@, suggesting a remarkable similarity in the 
processes governing star formation and dust production 
over the last 10 billion years of cosmic time. 

We have made significant progress in evaluating the 
effects of dust in typical star-forming galaxies at z ~ 2 
as traced by thermal infrared emission. However, even 
the superb sensitivity and resolution of Herschel are not 
sufficient to individually detect L* galaxies at these red- 
shifts. The most accurate estimate of dust luminosity, 
and hence total star formation rate, can only come from 
combining individual measurements of UV and IR lu- 



minosities. Though the observations will be targeted, 
ALMA promises to extend our knowledge of the dust 
continuum luminosities of individual star-forming galax- 
ies at z ~ 2, thus providing the key ingredient to com- 
bine with UV measurements. These developments point 
to a wealth of forthcoming information on the properties 
of dust in typical star-forming galaxies at high redshift. 
Ultimately, a robust characterization of the relationship 
between infrared and short wavelength (UV, Ha) emis- 
sion will be crucial for inferring dust attenuation and 
total star formation rates of the very faint and/or very 
high redshift galaxies that will be inaccessible to even 
the next generation of long wavelength observatories. 
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